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Abstract 

An asymptotic scaling theory is presented using the conceptual basis of trapping-free subspace 
(i.e., orthogonal subspace) to establish the generic mechanism of optimal efficiency of excitation 
energy transfer (EET) in light-harvesting systems. Analogous to Kramers' turnover in classical 
rate theory, the enhanced efficiency in the weak damping limit and the suppressed efficiency in the 
strong damping limit define two asymptotic scaling regimes, which are interpolated to predict the 
functional form of optimal efficiency of the trapping-free subspace. In the presence of static disorder, 
the scaling law of transfer time with respect to dephasing rate changes from linear to square root, 
suggesting a weaker dependence on the environment. Though formulated in the context of EET, the 
analysis and conclusions apply in general to open quantum processes, including electron transfer, 
fluorescence emission, and heat conduction. 
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The optimization of the excitation energy transfer (EET) process presents a challenge for 
understanding photosynthetic systems as well as for designing efficient solar energy devices. 
Multi-dimensional spectra have allowed detailed probes of EET dynamics, revealing signa- 



tures of quantum coherence 



-|4[]. In contrast to the common belief that noise retards motion 



and coherence enhances mobility, the EET efficiency can be optimized at an intermediate 



level of noise, leading to the notion of noise-enhanced energy transfer |5l4l0l|. Fermi's golden 
rule rate provides a simple interpretation and suggests that stochastic resonance between 
the donor and acceptor enhances EET efficiency |8l-ll0| . Anotherpossible mechanism is that 
noise suppresses destructive interference between pathways 9|. The situation is further 
complicated by the findings that initial preparation, coherence of incident photons, site en- 
ergy, spatial correlation, static disorder, and various approximations invoked in quantum 
master equations can all play a role in establishing optimal efficiency lQ-lsj]. Therefore, 
a general mechanism for optimization in an arbitrary EET system accompanying all these 
effects is clearly needed but has not yet been formulated. 

In this Letter, we utilize the concept of trapping-free subspace (i.e., orthogonal subspace) 
to bring together all of the above considerations into a unified framework, that allows us to 
establish asymptotic scaling relations under both dynamic and static disorder and identify 
the generic behavior of optimal EET efficiency. The definition of a trapping-free subspace 
is a generalization of the invariant subspace jg] and is closely related to the concept of 
decoherence-free subspace in quantum information 16 ]. 

Model. — We consider a light-harvesting EET system (see examples in Fig. [T]) described 
by the quantum master equation for the reduced density matrix of the single excitation 



manifold 



9] 



p(t) = -Cp{t). 



(1) 



Here, the Liouville superoperator C = C sys + £dissi P + Arap + £dcca y comprises four terms, 
each describing a distinct dynamic process: (i) £ sys p(t) = (i/h)[H$, p(t)], the dynamics 
of the isolated system, where H§ is the system Hamiltonian; (ii) £dissip> the exciton re- 
distribution and dephasing within the single-excitation manifold due to the interaction with 
the surrounding environment; (iii) £trap, the trapping of excitation energy at the reaction 
center for the production of chemical energy; (iv) £decay; the decay of the excitation energy 
to ground state in the form of heat or a photon. The Liouville superoperator £ becomes 
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invertible because of the two irreversible energy depletion terms, £ trap and £diss P - A basic 
property of an EET system is its energy transfer efficiency, 



/>oo 

q = Tr C tTauP p(t)dt, 
Jo 



(2) 



where Tr denotes the trace over states. For an efficient EET system such as a light-harvesting 
protein complex, its nearly unit energy transfer efficiency, q ~ 1, implies a clear time-scale 
separation between the decay process and the trapping process. Under this condition and 
with a homogeneous decay rate kd, the transfer efficiency is given, to a good approxima- 
tion by 

f 



q 



(3) 



where (t) = Tr [£" 1 p(0)] is the the average trapping time for the initial density matrix p(0) 
and Co = C sys + £ trap + £dissi P * s ^ ne new Liouville superoperator in the absence of decay. 
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FIG. 1: (color online) Various EET systems with the trapping free subspace: a) A donor-acceptor 
system with an energy mismatch A larger than the site-site coupling J (A > J), b) A homogeneous 
iV-site (N S> 1) chain where the site-site coupling decays exponentially with distance, c) A 
schematic diagram of a two-generation three- fold dendrimer [2] . 

The optimization of energy transfer efficiency is thus simplified to the minimization of 
(t). It has been found that for many EET systems, (t) decreases with increasing T for weak 
dissipation and diverges with T for strong dissipation, where T represents a characteristic 
dissipation strength in £diss P , e.g., the dephasing rate. [5r|l0|. A representative curve of this 



dependence is plotted in Fig. [2^, where the trapping time (t) is minimal at an optimal noise 
level r op t- To understand this generic behavior, we investigate the asymptotic scaling of the 
trapping time in the strong and weak dissipation limits. 

Asymptotics of the trapping time. - In the strong damping limit (T — > oo), quantum 
coherence is quickly destroyed by noise, and energy is transferred from donor sites to the 
trap through a series of incoherent hops. The hopping rate k^op can be estimated from 
Fermi's golden rule, k^op ~ l^| 2 /T, where J is the exciton coupling strength. This relation 
is consistent with the Forster theory where T is proportional to the spectral line-width. The 
average trapping time diverges linearly with T, giving 

(*>~*£i~r/|j| a . (4) 

In a recent paper, we mapped exciton dynamics to network kinetics using a systematic 
expansion and recovered hopping kinetics as the leading order term Q| in the strong damping 
limit. The exact T-dependence of the hopping rate fchop is usually more complicated for a 
quantum non-Markovian description of dissipation. Nevertheless, the linear T-dependence 
of the trap ping time in Eq. (J3J) can serve as a simple scaling relation for relative large values 



of damping [l7|] 



In the opposite weak damping limit (r — > 0), energy transfer can be enhanced by dynamic 
noise such that (t) decreases with increasing T. Here the starting point for studying EET 
dynamics is the delocalized exciton basis, i.e., eigenstates of Hs, and energy is transferred 
coherently to the trap state through delocalized excitons. However, not every exciton state 
is capable of efficient energy transfer; a set of excitons, orthogonal to the trapping operator, 
define the trapping- free subspace $j_ (i.e., orthogonal subspace), 

ArapP± = 0, (5) 

where the index _L denotes the trapping- free subspace such that p± G <E>j_. Without the 
coupling to the trap, Eq. ([5]) indicates that the population in $x is conserved and therefore 
the coherent trapping time diverges, < t > |r=p = oo. For a single trap site, the procedure 
of identifying $^ has been provided in Ref. 8[. We note that Eq. §5§ is analogous to the 
definition of decoherence-free subspace in quantum information, where the de-coupling from 



the system-bath interaction 16] is defined as the de-coupling from the system-trap overlap. 

The system-bath coupling induces interactions between the trapping-free and other exci- 
ton states, thus leading to population depletion from the trapping-free subspace. Since 



the population of the unperturbed orthogonal subspace will not deplete, dissipation of 
the orthogonal subspace $j_ dominates the population transfer time to the trap, giving 
(£o ~ (£dissi P ;_i_) _1 ) where C± represents the matrix element of the Liouville superoper- 
ator defined in For nonzero population in <E>j_, the leading order of the average trapping 
time is given by the survival time in the orthogonal exciton subspace, 

(t)^Tr[(/ dissip;X )^p±(0)]/r, (6) 

where the linear form £dissip ~ T^dissip is valid in the limit of T — > 0, and Idissp is the reduced 
Liouville superoperator, independent of V. 

Generality. - The two scaling relations in the asymptotic regimes are based on gen- 
eral physical arguments and thus independent of specific details such as the system-bath 
coupling, bath spectral density, truncation method, and approximations in the quantum 
master equation. Combining Equations (j4|) and (jHJ), we obtain an estimation of the optimal 
condition, 

ropt-iTr^^J-VxW]} 172 ^!, (7) 

which depends on the system Hamiltonian as well as the initial condition. In fact, the optimal 
efficiency or minimal trapping time is analogous to the Kramers turnover in reaction rate 
theory, where the two scaling regimes correspond to energy diffusion and spatial diffusion, 



respectively [20, \2l\ . The change of the reaction coordinate from energy to spatial position 
in classical rate theory corresponds to the change of basis set from excitons to local sites 
in energy transfer theory. However, the analogy to classical rate theory is limited to the 
orthogonal subspace and does not apply to the non-orthogonal subspace, where coherent 
energy transfer to the trap state is the dominant pathway in the weak damping limit and 
will not display noise-induced enhancement as illustrated in Fig. [2]A. 

The orthogonal condition in Eq. (jSJ) for $x can be realized in various systems, (i) In an 
inhomogeneous system with large energy mismatches (see Fig. la), the orthogonality can 
arise from a vanishingly small overlap coefficient between the donor and acceptor. This 
situation can be well described by Fermi's golden rule and qualitatively explains the optimal 



efficiency observed in FMO [lOl Il8|. (ii) In a spatially extended system with a large number 
of sites, the orthogonality can arise because the average overlap coefficient of a local site 
with the trap decreases with the system size, e.g., a long, homogeneous chain system (see 



Fig. lb), (iii) In a system with intrinsic symmetry, each exciton state resulting from the 
diagonalization of H sys is associated with a specific symmetry. A subset of these eigen-states 
are incompatible with the symmetry of the trapping Liouville operator £ tmp , thus leading 
to orthogonality. The fully-connected network 8| and the dendrimer in Fig. [Tb are examples 
of such topological symmetry. The orthogonality due to symmetry in case (iii) is rigorous, 
whereas the orthogonality in cases (i) and (ii) requires either a large energy mismatch (i.e. 
detuning) or a large system size. In cases (i) and (ii), the average trapping time may not 
diverge at T = but can still lead to an optimal efficiency because of noise-enhancement 
under weak non-orthogonality. 
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FIG. 2: EET in the dendrimer depicted in Fig. a) The trapping time (t) vs. the pure dephasing 
rate T. For the incoherent initial condition defined at six outer 2nd-generation sites (i.e, the non- 
trapping state), the solid line is a full calculation. The dotted-dashed line is the large-r result from 
Eq. ([3]). The dashed line is the small-r result from Eq. ©. As a comparison, the dotted line is 
the result for the coherent initial condition defined at three lst-generation sites (i.e, the trapping 
state), b) The transfer efficiency q vs. T for the 1st initial condition. The dashed line is the exact 
result as calculated from Eq. (T2]) while the solid line is obtained from the approximation in Eq. ([3]) 
and modification for T < 1. Here parameters are J = 20 meV, kt = 5 meV, and kd = 5 /xeV. 
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An example — To verify the above analysis, we consider a two-generation three-fold den- 
drimer depicted in Fig. jl9]. A tight-bonding model is used for the system Hamiltonian, 
giving Hs-ij = EiSij\i)(i\ + — Sij)\i){j\, where \i) represents a localized excited state at 
site i. All the site energies e$ are assumed to be identical and the site-site interaction = 20 
meV is defined between two connected sites. The center site is the trap where the trapping 
process occurs with rate kt = 5 meV. The decay process is characterized by a homogeneous 
rate kd = 5 /ieV. For simplicity, we approximate dissipation by the Haken-Strobl-Reineker 



(HSR) model [22|, |23( and consider homogeneous pure dephasing, £dissip;y = (1 _ 
where T is defined as the pure dephasing rate for each coherence element pij. For this den- 
drimer system, the symmetric Hamiltonian has an orthogonal subspace consisting of seven 
exciton states. We now compare two different cases of initial preparation. In the first case, 
an incoherent population is evenly distributed at six outer 2nd-generation sites, leading to 
P_l(0) 7^ 0. The average trapping time is plotted as a function of the pure dephasing rate 
in Fig. PSti. The divergence of (t) in the strong and weak dephasing limits is in an excellent 
agreement with the asymptotic behaviors predicted in Eqs. (OS) and ([6]), respectively. Fig- 
ure [2b shows that our approximate equation in Eq. (J3j) provides a quantitatively accurate 
description for the EET efficiency. Here a slight modification is applied in Eq. ([3]) to correct 
the zero-dissipation transfer efficiency 2Jj. In the second case, a coherent state is evenly 
distributed at three lst-generation sites, leading to p_i_(0) = 0. The small-r divergence of (t) 
in Eq. ([6]) disappears since no initial population exists in the orthogonal subspace. Above 
a threshold at the intermediate dephasing rate, (t) changes from a plateau to the same lin- 
early increasing function of T. This calculation confirms that dynamic noise can enhance the 
EET if the exciton subspace $j_ is orthogonal to the trapping process, and demonstrates the 
relevance of the initial condition for the energy transfer efficiency and thus the possibility of 
quantum control of energy transfer efficiency. 

Static disorder. - - We introduce an energy disorder 5ei at each site i, that follows the 
Gaussian distribution, P{5ei) = exp[— (5ef /(2af))]/ \f 7 hxo i , where Oi is the variance of energy 
disorder. The resulting ensemble average is given by (x) a = Hi J x(5ei)P(5ei)d5ei with 
x = (t) or q. Using the above dendrimer model with the first initial condition, we present the 
results of {{t)) a and (q) a obtained from a Monte Carlo simulation of 10 5 samples. As shown in 
Fig.[3li, static disorder is irrelevant in the strong dephasing limit (r 3> cr). A more interesting 
a-dependence occurs in the weak dephasing limit (r a), where the trapping time is 
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a (meV) T (meV) 

FIG. 3: The ensemble-averaged EET for the same dendrimer as in Fig. [2] with a static disorder 
a = 4 meV. a) The solid line is the simulation result of {{t)) a , referenced by the dotted line without 
a from Fig. [2^. The dashed line is the dominant contribution from the orthogonal subspace, and 
the dashed-dotted line is the fitting result of 1/VT. b) {q) a vs. T: the solid line is the ensemble 
average of Eq. ([2]) whereas the dashed line is from Eq. ([3]) using {{t)) a and modification for T < 1. 
c) The weak-dephasing u-dependence of ((t)) a , with F = 10~ 12 , 1CT 11 , 10~ 10 , 1CT 9 , and 1CT 8 meV 
(from top to bottom). Symbols denote simulation results, whereas the solid lines are fitted with 
{{t))a = c' /a. d) The fitting coefficient c' (circles) is plotted as a function of T and can be further 
fitted by d = c/VT (the solid line). 
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dominated by the survival time in the trapping-free subspace. Static disorder can destroy the 
orthogonality in Eq. (jSJ) and induce a large reduction of the trapping time. However, if the 
random energy disorder falls below the amplitude of stochastic noise, the weak orthogonality 
is preserved and the trapping time diverges as (t) oc 1/T. Since the amplitude of the 



where Se describes the effective in-phase energy fluctuation over all the sites, P(5e) ~ \ ja is 
the probability distribution, and (t)\\s e \ <a > ~ (t)\s £ =o ~ 1/r is approximately uniform within 
this regime. The pre-factor c depends on the initial condition, the system Hamiltonian, and 
the trapping rate. We confirm the asymptotic relation, ({t)) a oc 1/y/T, by numerical fitting 
in Fig. [3K Further, we calculate the weak-dephasing ((t)) a as a function of a for different 
values of T in Figs. [3b-d and rigorously estabalish the scaling relation predicted by Eq. (jSj), 
i.e., ({t)) a oc 1/cr. Our calculations suggest that nature can use both static and dynamic 
disorders cooperatively to achieve efficient and robust energy transfer. Indeed, as shown in 
Fig. [3b, the zero dephasing efficiency (g|r=o)<r is drastically enhanced from 0.2 to 0.8 with 
a = J/5. 

Conclusion. - In this Letter, we demonstrate that the generic mechanism of noise 
enhanced EET is to assist energy flow out of the orthogonal exciton subspace, and the 
competition between noise-enhanced EET in the weak dissipation regime and noise-induced 
localization in the strong dissipation regime leads to an optimal efficiency. We determine the 
scaling relations of the average trapping time in these two regimes and use the asymptotic 
relations to qualitatively predict the optimal noise. The presence of static disorder reduces 
the exponent of divergence in the weak-dissipation limit and thus makes the EET process 
more robust against noise. Our analysis is not limited to EET but also applies in general to 
electron transfer, fluorescence emission, heat conduction and other open quantum processes. 
Sophisticated numerical methods, such as the hierarchy equation approach, will be used to 
quantify the general optimal condition in the EET process. 
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stochastic noise is proportional to a/T, the range of relevant random energy that contributes 
to the divergence is defined by o' ~ vT. Thus, we obtain the new asymptotic relation in 
the weak dephasing limit by an integral over these small disorders, giving 
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